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Abstract 

Traditional treatments of the first law of black hole thermodynamics do not 
include a discussion of pressure and volume. We give an overview of recent 
developments proposing a definition of volume that can be used to extend 
the first law to include these appropriately. New results are also presented 
relating to the critical point and the associated second order phase transition 
for a rotating black-hole in four-dimensional space-time which is asymptot- 
ically anti-de Sitter. In line with known results for a non-rotating charged 
black-hole, this phase transition is shown to be of Van der Waals type with 
mean field exponents. 
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1 Introduction 



Ever since Hawking's discovery in 1974, [TJ |2j [3], that black holes have a 
temperature associated to them, in the simplest case a temperature inversely 
proportional to their mass, 

(we use units in which c = 1), the thermodynamics of black holes has been 
a fascinating area of research. Equation (JT|) immediately implies that a 
Schwarzschild black hole in isolation is unstable: it will radiate and in so 
doing loses energy hence the mass decreases, thus increasing the tempera- 
ture causing it to radiate with more power leading to a runaway effect. 

Hawking's result is fundamentally quantum mechanical in nature and 
came after a number of important developments in the classical thermody- 
namics of black holes. Penrose [3] realised that the mass of a rotating black 
hole can decrease, when rotational energy is extracted, and this was followed 
by the observation that the area never decreases in any classical process. 
Nevertheless there is still a minimum, irreducible, mass below which one 
cannot go classically |6|. This lead Bekenstein's to propose that an en- 
tropy should be associated with a black hole that is proportional to the area, 
A, of the event horizon, [S] and Hawkings result ([T]) fixes the co-efficient 
to be one-quarter. In natural units 

where TiG is the Planck length squared. 

The first law of black hole thermodynamics, in its simplest form, asso- 
ciates the internal energy of a black hole with the mass, U (S) = M, (more 
precisely the ADM mass, as defined with reference to the time-like Killing 
vector at infinity [9]) and reads 

dU = TdS. (3) 

The black hole instability referred to above is reflected in the thermody- 
namic potentials by the fact that the heat capacity of a Schwarzschild black 
hole, 

C-T^---^<0 (A) 

is negative. 
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The first law generalises to electrically charged, rotating black holes as 

dU = TdS + VLdJ + <$>dQ (5) 

where J is the angular momentum of the black hole, Q its angular velocity, 
and Q the electric charge and the electrostatic potential (see e.g. [9]). 

In contrast to elementary treatments of the first law of black hole ther- 
modynamics it is noteworthy that (jSJ) lacks the familiar PdV term and there 
has recently been a number of papers addressing this issue, [TOj [HI [12j [13j EHJ 
IT5| [T6| ITT] . This paper gives an overview of these developments and presents 
some new results in this regard, relating to the presence of a critical point and 
a second order phase transition for a rotating black-hole in four-dimensional 
asymptotically anti-de Sitter space-time, a critical exponent that is of a Van 
der Waals type with mean field exponents. 

A little thought shows that it is by no means obvious how to define the 
volume of a black hole. For a Schwarzschild black hole the radial co-ordinate, 
r, is time-like inside the event horizon, where r < r^, so it would seem non- 
sensical to associate a volume V = 4n J^ h r 2 dr = with the black hole. 
In fact identifying any function of alone with a volume, V{rh), will lead 
to inconsistencies in a thermodynamic description since the area, and hence 
the entropy, is already a function of r^, S = nrl, so any volume V{jh) would 
be determined purely in terms of the entropy. The internal energy, U(S, V), 
should be a function of two variables, so giving V(S) uniquely as a specific 
function of S is liable to lead to inconsistencies. We shall see below how this 
potential problem is avoided. 

2 Pressure and Enthalpy 

From the point of view of Eintein's equations a pressure is associated with 
a cosmological constant. There is now very strong evidence that the cosmo- 
logical constant in our Universe is positive [18j [19]. This poses a problem 
for the study of black hole thermodynamics for two reasons: firstly there 
is no asymptotic regime in de Sitter space which allows the unambiguous 
identification of the ADM mass of a black hole embedded in a space with a 
positive A; secondly positive A corresponds to negative pressure, implying 
thermodynamic instability. The first problem is related to the fact that there 
are two event horizons for a de Sitter black hole, a black hole horizon and a 
cosmological horizon, and the radial co-ordinate is time-like for large enough 



2 



values of r, outside the cosmological horizon. The second problem is not 
necessarily too serious as one can still glean some information from nega- 
tive pressure systems which are thermodynamically unstable [20J (instability 
is not an insurmountable barrier to obtaining physical information from a 
thermodynamic system, after all, as described above, Hawking's formula (PQ) 
shows that black holes can have negative heat capacity but it is still a central 
formula in the understanding of black hole thermodynamics). In contrast for 
negative A there is no cosmological horizon and the pressure is positive, the 
thermodynamics is perfectly well defined, so we shall restrict our considera- 
tions here to negative A and identify the thermodynamic pressure P = —-^q 
with the fluid dynamical pressure appearing in Einstein's equations. 

The notion that the cosmological constant should be thought of as a 
thermodynamic variable is not new, and its thermodynamic conjugate is 
often denoted 9 in the literature, [2U [22l [231 EH EH [26l [271 12H1 [291 1301 , but 
was not given a physical interpretation in these works. 

It may seem a little surprising to elevate A to the status of a thermody- 
namic variable. A is usually thought of as a coupling constant in the Einstein 
action, on the same footing as Newton's constant, and it might seem bizarre 
to think of Newton's constant as a thermodynamic variable. However the 
nature of A has long been mysterious [31] and we should keep an open mind 
as to its physical interpretation. Indeed in [TT] it was argued that A must be 
included in the pantheon of thermodynamic variables for consistency with 
the Smarr relation |32j, which is essentially dimensional analysis applied to 
thermodynamic functions. Furthermore [11] suggested that, for a black hole 
embedded in anti-de Sitter (AdS) space-time, the black hole mass is more 
correctly interpreted as the enthalpy, H beloved of chemists, rather than the 
more traditional internal energy, 

M = H(S, P) = U(S, V) + PV. (6) 

The PV term in this equation can be though of as the contribution to the 
mass-energy of the black hole due the negative energy density of the vacuum, 
e = — P, associated with a negative cosmological constant. If the black hole 
has volume V then it contains energy eV = —PV and so the total energy is 
U = M -PV. 

This interpretation forces us to face up to the definition of the black 
hole volume. In [TT] V is defined as the volume relative to that of empty 
AdS space-time: the black hole volume is the volume that is excluded from 
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empty AdS when the black hole is introduced. We shall refer to this as the 
"geometric volume" below. Other suggestions for the volume of a black hole 
have been made in [33], [31] 

An alternative definition of the black-volume is that it is the thermody- 
namic conjugate of the pressure, under the Legendre transform (jSJ), 

which we shall call the "thermodynamic volume". 

With the definition of the thermodynamic volume (j7|) we are in a position 
to state the definitive version of the first law of black hole thermodynamics, 

dU = TdS + VtdJ + <$>dQ - PdV (8) 

which follows from the Legendre transform of 

dM = dH = TdS + VLdJ + $rfQ + VdP. (9) 

Equation ([9]), in dA notation, appeared in [35]. A PdV term was considered 
in [36] and [12], but only for J = 0, when V and S cannot be considered to 
be independent thermodynamic variables. 



3 Thermodynamic volume 

The suggested definition of the thermodynamic volume ([7j) must be tested for 
consistency. For example, for a non-rotating black hole in four-dimensional 
space-time, the line element is given, in Schwarzschild co-ordinates, by, 

d 2 s = -f(r)dt 2 + r\r)dr 2 + rW, (10) 

with 

„/ \ 2m A 9 . , 

f(r =1 ^r 2 } (11) 

r 3 

and dVt 2 = d6 2 + sin 2 9d(f> 2 the solid angle area element Q The event horizon 
is defined by /(r^) = 0, 

|rg-r h + 2m = 0, (12) 

^rom now on we set G = fi = 1 to avoid cluttering formulae. 
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but we do not need to solve this equation explicitly in order to analyse (|7J). 
We already know that 



S = n$, P = — (13) 

and, for negative A, the ADM mass is M = m [23], which, following the 
philosophy of [IT], we identify with the enthalpy, H(S,P). Solving ( 112]) for 
m immediately yields 



m = f I 1 - rv ' (14) 

from which H(S, P) — M — m, with f[T3"j) . identifies the enthalpy as 

The usual thermodynamic relations can now be used to determine the 
temperature and the volume, 

That the resulting thermodynamic volume (for a non-rotating black hole) 
is identical to the geometric volume is quite remarkable, but appears co- 
incidental as this equality no longer holds for rotating (Kerr-AdS) black holes, 
as we shall see. It does however hold for non-rotating black holes in all 
dimensions [36], [121 H3] • 

As mentioned in the introduction, equation (|T7|) has a potential problem 
associated with it, in that it implies that the volume and the entropy cannot 
be considered to be independent thermodynamic variables, S determines V 
uniquely - they cannot be varied independently and so V seems redundant. 
Indeed this may the reason why V was never considered in the early literature 
on black hole thermodynamics. But this is an artifact of the non-rotating 
approximation, V and S can, and should, be considered to be independent 
variables for a rotating black hole. 
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The line element for a charged rotating black hole in 4- dimensional AdS 
space is [37] 

A / asin 2 6 , \ 2 p 2 j2 p 2 - A,sin 2 #/ , r 2 + a 2 \ 2 

rfs 2 = — - dt — deb +^-dr 2 +^—d6 2 +^^: adt — dd> 

P 2 \ = / A A e P 2 \ - J 

(18) 

where 

(r 2 + a 2 )(L 2 + r 2 ) 2 a 2 , „ 

A = ^ ^ ^-2mr + g 2 , A„ = 1 - — cos 2 #, 

p 2 = r 2 + a 2 cos 2 e, S = l-^, (19) 

and the cosmological constant is A = — -p- = — 87rP. 

The physical properties of this space-time are well known |23j. The metric 
parameters m and q are related to the ADM mass M and the electric charge 
Qby 

M=^, Q = |. (20) 

The event horizon, r + , lies at the largest root of A(r) =0, so, in terms of 
geometrical parameters, 

M _ (r 2 + + a 2 )(L 2 + r 2 + ) + q 2 L 2 

2r,L 2 E 2 1 ' 



r 2 + a 2 



and the area of the event horizon is 

A = 4tt ' + _ - , (22) 
giving 

S = *Z±±. (23) 

The angular momentum is J = aM and the relevant thermodynamic 
angular velocity is 

n a < L2 + '"> (24) 

As explained in [35J , Q here is the difference between the asymptotic angular 
velocity and the angular velocity at the black hole outer horizon. 
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The electrostatic potential, again the difference between the potential at 
infinity and at the horizon, is 



* = (25) 
r± + a 2 



To determine the thermodynamic properties, M must be expressed in 
terms of S, J, Q and P (or, equivalently, L). This was done in |35j and the 
result is 



H(S,P,J,Q) :=-\ 



(26) 



This generalises the Christodoulou-Ruflini formula [6] for the mass of a 
rotating black hole in terms of its irreducible mass, Mi„. (The irreducible 
mass for a black hole with entropy S is the mass of a Schwarzschild black 
hole with the same entropy, M 2 rr — j-)- 
The temperature follows from 



dH 
~dS 

1 

8nH 



(27) 



J,Q,P 
1 



tiQ 2 8PS" 



1 



5 



5 3 

from which we immediately see that T > requires 
S 2 (_ nQ 2 8PS\ ( _ nQ 2 



8PS 



^ 1 j 



r < 



An 2 



nQ 2 8PS" 

1 + — + — 



1 



S 



8PS 



(28) 



The maximum angular momentum, 

J max 



S_ 

2n\ 



tiQ 2 8PS 

1 + -T7- + —=~ 



vrQ 2 



s 



s 



+ 8PS 



(29) 



is associated with an extremal black hole. 

From (J7|) and fl26l) the thermodynamic volume is 



V 



dH 



S,JA 



3nH 



S [S + ttQ 2 + 



8PS 2 



+ 2n 2 J 2 



(30) 
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which is manifestly positive. 

The angular velocity and the electric potential also follow from 



n 



dH 
~dJ 



S,Q,P 



4vr V(l + ^ 



and 



OH 
dQ 



2nQ (S + ttQ 2 



8PS 2 



S,J,P 



2HVnS 



via 



(31) 



(32) 



The Smarr relation follows from (gBJ), (1271), (13Uj), (13Tj) and 



namely 



H 



+ PV-ST-JQ 



0. 



(33) 



from which it is clear that the PV-term must be included for consistency, as 
pointed out in [TT| . 

It is clear from fl3"U|) that, in general, V is a function of all the four 
independent thermodynamical variables, S, P, J and Q, but for the limiting 
case J = 0, 



V 



4Sl_ 
3 a/?' 



(34) 



is determined purely in terms of S alone, independent of both P and Q. 
Thus, as explained in the introduction, V and S cannot be viewed as ther- 
modynamically independent variables as J — > 0, rendering the description in 
terms of the thermodynamic potential U(S, J) impossible in this limit. 

Expressing the thermodynamic volume (1301) in terms of geometrical vari- 
ables one gets [15] 



V 



2vr j (r 2 + + a 2 )(2r 2 + L 2 + a 2 L 2 - r 2 + a 2 ) + L 2 q 2 a 2 



Given that the area of the event horizon is 

A = Att 

then, if we define a naive volume 

Vo 



ri + a 2 



r+A 47r r + (r 2 + a 2 ) 



3 
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(35) 



(36) 



(37) 



equations (|2"T|) and (ESJ) give 



V=V + 



4vra 2 M 




(38) 
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a formula first derived in [14J. As pointed out in that reference, equation 
fl37|) implies that the surface to volume ratio of a black hole is always less 
than that of a sphere with radius r + in Euclidean geometry. This is the op- 
posite of our usual intuition that a sphere has the smallest surface to volume 
ratio of any closed surface — the isoperimetric inequality of Euclidean ge- 
ometry. Thus the surface to volume ratio of a black hole satisfies the reverse 
of the usual isoperimetric inequality (a similar result holds in higher dimen- 
sions [H]). At least this seems to be the case if quantum gravity effects are 
not taken into account. In one case where quantum gravity corrections can 
be calculated using the techniques in [38J , the three-dimensional Banados- 
Zanelli-Teitelboim (BTZ) black hole [39], they tend to reduce the black hole 
volume [13J so it seems possible that quantum gravity effects may affect the 
reverse isoperimetric inequality. 



To examine the consequences of the PdV term in the first law we need to 
perform a Legendre transform on the enthalpy to obtain the internal energy 
U(S, V, J, Q) from U = H — PV. We first write the enthalpy (126]) in the form 



4 The First Law 



H= Va + bP + cP 2 , 



(39) 



where 



a 



b 



c 




(40) 



Note that the discriminant, 




(41) 
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is positive. 
Now 



V 



OH 



S,J,( 



b + 2cP 
2H 



2HV -b 
2c ' 



This allows us to re-express H function of V, 



H 



1 b 2 - Aac 



2 V V 2 -c 



We can immediately conclude that 



with equality only when 



b - Aac 



64vr 2 



-J 2 J 2 



9 



SQ 2 



71 



i.e. when J = 0. 

It is now straightforward to determine 



U = H -PV = H 



HV 2 bV\ bV 



{V 2 - c){b 2 - Aac) 



2c / 2c 



2c 



which immediately gives 



U(S,V,J,Q) 



3V 
An 



~\J\ 



S 
2n~ 



S 



7T 



+ Q ) + J 



~An~) 



71 



31 \ 



Note the subtlety in the J — y limit, ( I47jl is not differentiable at J 
unless 



/3^x2 _ ,S^3 



there. 
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Equation (H7|) can now be used to study the efficiency of a Penrose process. 
If a black hole has initial mass M iy with internal energy Ui, and is taken 
through a quasi-static series of thermodynamic steps to a state with final 
internal energy Uf, then energy can be extracted if Uf < Ui. This is the 
thermodynamic description of a Penrose process [I] and the efficiency is 

Ui -u f 

" = ~M~' (49) 

We can determine the maximum efficiency for a process at constant P by 
first expressing U in (137|) in terms of 5, P, J and Q: 

(5 + ttQ 2 ) (5 + ttQ 2 + + 47T 2 (l + U£) J 2 

U = V j M 37 V 3 i^, (50) 



5 + ^2+8^)^4,2 (l + ^) P 



which is manifestly positive. 

For simplicity consider first the Q = case, for which 

dU = TdS + VldJ - PdV. (51) 

The work extracted at any infinitesimal step is 

dW = -dU = -TdS - VLdJ + PdV (52) 

and, since dS > 0, this is maximised in an isentropic process dS = 0. Now 
with Q = and S and P held constant, the internal energy in equation 
f )50|) can be thought of as a function of J only, U(J). The greatest efficiency 
is then obtained by starting with an extremal black hole and reducing the 
angular momentum from J m ax to zero, it is given by 

u(J max ) - u(Q) 

Vext — T-7- r {06) 

{J max ) 

where H(J max ) = M ext is the initial extremal mass. One finds 

1 + 2PS _ __L 3_ 

Vext 1 + Aps ^2 + 8PS{3 + 8PS)' [ ' 

In asymptotically flat space, A = 0, we set P = in r] ext and obtain the 
famous result [9] 

Vext = 1 - i w 0.2929. (55) 
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More generally, r] ext is a maximum for SP = 1.837 . . . (obtained by solving a 
quartic equation) and attains there the value 0.5184 .... Thus turning on a 
negative cosmological constant increases the efficiency of a Penrose process, 
as first observed in [15] . 

What is happening here is that, as \J\ decreases (giving a positive con- 
tribution to dW) the volume decreases, which actually tends to decrease the 
work done because of the PdV term in (15121 . But when P > 0, the extremal 
value | Jmax\ in (129]) is increased, which more than compensates, and overall 
rj ext is increased. 

For a charged black hole the internal energy is a function of J and Q for 
an isobaric isentropic process, U(J,Q). The requirement J max > in (1291) 
imposes the constraint 

Q 2 <QL x = (^)a + 8PS) (56) 

on the charge. The greatest efficiency is achieved starting from an extremal 
black hole with Q 2 = Q 2 max and reducing both J and Q to zero in the final 

U( Jmax, Qmax) ~ J(0, 0) _ 3 / 1 + 8PS \ 
H^Jmaxi Qmax 

) 2V3 + 16PS/' v ' 

with H(J max , Qmax) the initial extremal mass, M ext . For large S efficiencies 
of up to 75% are possible [15], which should be compared to 50% in the A = 
case, [9]. 



5 Critical Behaviour 

With knowledge of both H and U general questions concerning the heat 
capacity of black holes can be addressed. The heat capacity at constant 
volume, Cy = Tj (f|) VJ > tends to zero when J = 0, though Cy can be 
non-zero for J ^ it does not diverge. For comparison the heat capacity at 
constant pressure, Cp = T / j pj , Cp vanishes when T = and diverges 

when |^ =0. 

A full stability analysis was given in |35j and there are both local and 
global phase transitions. Local stability can be explored visually, by plotting 
thermodynamic functions, or analytically, examining the curvature of the 
derivatives of thermodynamics functions. 
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5.1 Q = 



Let us first focus on the Q = case. The red (lower) curve in the figure 
below shows the locus of points where Cp diverges in the J — S plane, it is 
given by setting the denominator of Cp, 

144 (tt JP) 4 (9 + 32 SP) + 24 (ttPJ) 2 (PS) 2 (3 + 16 SP) (3 + 8 SP) 2 

-(PS) 4 (1 -8 SP) (3 + 8 SP) 3 , (58) 

to zero. The green (upper) curve is the T = locus, all points above and 
left of this curve are unphysical as T < in this region. 



0.01 
0.008 
0.006 

JP 

0.004 



0.002 

o r 

0.02 0.04 0.06 0.08 0.1 0.12 0.14 

SP 

Figure 1: T = and Cp — > oo curves in J — S plane. 

There is also a global phase transition, not shown in the figure, when the 
free energy of pure AdS is lower than that of a black hole in asymptotically 
AdS space-time, the famous Hawking-Page phase transition |40]. We shall 
focus on the second order local phase transition here and examine its critical 
properties. 

In general, at fixed P and J, there are two values of S at which Cp 
diverges, and there is a critical point where these two values coalesce into 
one, the maximum of the lower curve in figure 1. This critical point was first 
identified in [35]. On purely dimensional grounds PCp can be expressed as 
a function of PS and PJ and the critical point can be found analytically, 
by solving a cubic equation, but the explicit form is not very illuminating. 
Numerically it lies at 
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« 0.7811 (60) 

crit 



(PS) crit w 0.08204, (PJ)crit ~ 0.002857. (59) 
The critical temperature is obtained from (j27j) . with Q — 0, 

/ T ' 
T P 

and the critical volume likewise from (130]) 

(^ 3/2 ) cr ,. t ~ 0.01768 (61) 

(the authors of [35J fix P = ^ « 0.1194, corresponding to L = 1, and find a 
critical value of J at J c w 0.0236). 

The equation of state cannot be obtained analytically, but its properties 
near the critical point can be explored by a series expansion and critical 
exponents extracted. Define the reduced temperature and volume as 

T - T c V -V c 

t = v = — -A (62) 

It is convenient to expand around the critical point using 

p:= 16n(PJ- (P J) crit) (63) 



and 

q:=8(PS-(PS) cnt ). (64) 

Expanding the temperature (!27l) around the critical point, with Q set to zero, 
gives 

t = 2.881 p + 2.201 pq + 0.3436 q 3 + o(p 2 ,pq 2 , q 4 ). (65) 
while similar expansion of the thermodynamic volume (130 j) yields 

v = -10A4p + 2.28Aq + o{p 2 ,pq,q 2 ). (66) 

For a given fixed J > 0, p is the deviation from critical pressure in units of 
l/(167r J), but one must be aware that this interpretation precludes taking 
the J — > limit in this formulation. Bearing this in mind, f )65|) and fl66|) give 
the J > 0, Q = equation of state parametrically in terms of q. Eliminating 
q one arrives at 

p = 0.3472 1 - 0.1161 tv - 0.02883 v 3 + o(t 2 , tv 2 , v 4 ). (67) 



14 



The critical exponent a is denned by 

c v oc r a 



(68) 



and, since as already stated, Cy does not diverge at t — 0, a — 0. To see 



this explicitly note that Cy = T/^\ v 



and 



dT 
dS 



v 



dT 
dS 



+ 



dT 
dP 



s 



dP 

as 



V 



c \dS 



+ 



dt_ 

dP 



OP 

as 



(69) 



Now, near the critical point, (1651) gives 

= sp dt 

I' 



dS 

dt_ 
dP 



dq 



o(p, q 2 



8 4 

oq 



op 



2.881(167r J) +o{p,q), 



(70) 
(71) 



implies dp = 0.2188 dq for constant v, from which is follows that 



while 

dp 

dS 

critical point and so a = 0. 

The exponent (3 is defined by 



is non-zero at the critical point, hence |^ 



does not vanish at the 



(72) 



where v > is the greater volume and v < the lesser volume across the phase 
transition, at constant pressure, when t < (v < is negative, since t> = at 
the critical point). 

Keeping p and t constant in fl67|) implies that 



P 



f V> dv = 0.3742 1 f V> dv- H (o. 1161 tt» + 0.02883 v 3 W (73) 



Allowing for the area of the rectangle in figure 2, namely 0.3742 |t|(i>> — f<), 
Maxwell's equal area law then requires 



v> 



(0.1161 to + 0.02883 v 3 )dv = 



Itlcxivl+vl). (74) 



It is clear from the figure that v > — v < ^> v > + f< so 



(75) 
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p 




Figure 2: Construction associated with Maxwell's equal area law. 
giving 

\t\ oc 0> - v<) 2 (76) 

and — ~. 

The critical exponent 7 is related to the isothermal compressibility, 

1 (dv\ _ 1 (dv\ (dv\ {wis j 
= ~v \op) T j = -v \ep) SJ - Kes) PJ (f)^ (77) 

which diverges along the same curve as Cp does (the adiabatic compression, 
k s = — y ywjgj' ^ S ever J rwnere finite — see equation fl89l) ). 7 gives the 
divergence of the isothermal compressibility near the critical point, 

k t oc r 7 . (78) 

7 can be found by expanding the denominator of Cp in (|58|) around the 
critical point, but a quicker method, since we have the equation of state, is 
to differentiate (1671) with respect to v, keeping t constant, giving 



dp 

dv 



oc -t, (79) 
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hence 



Kt oc 



dv 
dp 



1 



(80) 



and 7 = 1. 

Lastly setting t = in fl67|) we see that 

|p| °c M* 5 



(81) 



with 5 = 3, again the mean field result. 
To summarise, the critical exponents are 




(82) 



These are the same critical exponents as the Van der Waals fluid and, more 
importantly, are mean field exponentsJl 

It was first pointed in [321 [32] that a non-rotating, charged black hole has 
a critical point of the same nature as that of of a Van der Waals fluid, and 
the critical exponents for the black hole phase transition in this case were 
calculated in [T7] and verified to be mean field exponents, which are indeed 
the those of a Van der Waals fluid. A similarity between the neutral rotating 
black hole and the Van der Waals phase transition was first pointed out in 
[35] and further explored in [15] . 

The critical point can be visualised by plotting the Gibbs free energy 



for J = 1 and Q = 0, as a function of P and T as in figure 3. We see 
the "swallow-tail catastrophe" that is typical of the Van der Waals phase 
transition [33]. 

This structure is a straightforward consequence of Landau theory, |45j. 
Near the critical point the Landau free energy is 



L(T,P,v) = G(T,P) + A{(p-Bt)v + Ctv 2 + Dv 4 } + ... , (84) 



where G(T, P) is the Gibbs free energy and A, B, C and D are positive 
constants (for simplicity the constant J is not made explicit). As stressed 

2 While this manuscript was in preparation we became aware that these critical expo- 
nents have also been found, using a virial expansion approximation, in [41 . 



G(T, P, J) = H(S, P, J) - TS, 



(83) 
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Figure 3: Gibbs free energy as a function of pressure and temperature, at 
fixed angular momentum. 

in [36] L is not strictly speaking a thermodynamic function as it depends on 
three variables, p, t and v instead of two: v is to be determined in terms of 
p and t by extremising L to obtain the equation of state. 

For notational convenience equation fl84l) can be written, for fixed p and 
t, as 

L = a + bv + cv 2 + v 4 (85) 

where a, b and c need not be positive and L — > -r^L has been rescaled by a 
trivial positive constant. We are to think of b and c are control parameters 
that can be varied by varying p and t. 

Extremising (|85|) with respect to v determines the value of v in terms of 
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b and c through 



b = -lev - 4v 3 . 



(86) 



Using this in L leads to 

L = a-cv 2 - 3v 4 . (87) 

Equations fl86|) and ( IHTj) together give L(a, b, c) implicitly: a parametric plot 
of L(b, c), for any fixed a, reveals a characteristic "swallow-tail catastrophe" 
structure. With hindsight the swallow-tail structure is clear: in the A — D — 
E classification of critical points of functions, [J7], ( 185|) has three control 
parameters and is derived from type A\ in Arnold's classification. 

5.2 Q^O 

The above structure was first found in AdS black hole thermodynamics in the 
charged J = case [12J S3], where the equation of state can be found exactly 
and the critical exponents can be determined [17]. When both J and Q 
are non-zero an analytic analysis is much more difficult, for example finding 
the zero locus of the denominator of Cp requires solving a quintic equation. 
However numerical studies show that for a charged rotating black hole, as 
long as the charge is below the extremal value, the picture is qualitatively 
the same: the critical exponents are the same, the Landau free energy is still 
related to type A 4 and the Gibbs free energy still takes on a characteristic 
swallow-tail shape. For fixed values of J and Q, not both zero, all that 
changes is the numerical value of the co-efficients in equations (16"5"j) , (ESJ) and 
f )67j) or, equivalently the numerical values of the constants A, B, C and D in 
As long as none of these constants actually changes sign the nature of 
the critical point does not change and the critical exponents are the same. 

As first observed in [32] there is a line of second order critical points in 
the J — Q plane, as shown below, 

As long as J and Q are not both zero conservation of charge/angular mo- 
mentum protects the black hole against the Hawking-Page phase transition 
to pure anti-de sitter. 

6 Compressibility and the speed of sound 

In the previous section, the nature of the singularity in the isothermal com- 
pressibility near the critical point was discussed, but the adiabatic compress- 
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Figure 4: The line of second order critical points in the J — Q plane, [35] 



ibility 



V \dP 



(88) 



T,J,( 



is also of interest, and this was studied in [16] on which most of this section 
is based. From (J7J) one finds, setting Q = for simplicity, that 



36(2vrJ) 4 5 



(3 + 8PS){(3 + 8PS)S 2 + (2vr J) 2 }{2(3 + 8PS)S 2 + 3(2vr J) 2 } ' 



(89) 



This is finite at the critical point, indeed it never diverges for any finite 
values of S, P and J, and it vanishes as J — > 0: non-rotating black holes are 
completely incompressible. Black holes are maximally compressible in the 
extremal case T = 0, when J = J max in ( 1291) . 



2S{1 + 8PSf 



(90) 



T "° (3 + 8PS) 2 (l +APS)' 
A speed of sound, cs, can also be associated with the black hole, in the 
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usual thermodynamic sense that 
-2 d P 



dP 



9 2vrJ 4 , . 

= l+ p Ks = l + i >- _ 91 

SJ H {2(3 + 8PS)S 2 + 3 (2ttJ) 2 } 2 V ; 



where p = y is the density, cs is unity for incompressible non-rotating black 
holes and is lowest for extremal black holes in which case 



c s 



t=o \3 + 8PS 



giving c| = 0.9 (in units with c = 1) when P = 0. In the limiting case 
PS — > oo, Cg achieves a minimum value of 1/2. 

These results show that the equation of state is very stiff for adiabatic 
variations of non-rotating black holes and gets softer as J increases. For 
comparison, the adiabatic compressibility of a degenerate gas of N relativistic 
neutrons in a volume V at zero temperature follows from the degeneracy 
pressure 
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p -=^t(W ^ K °=4k, (93) 

For a neutron star — w 10 45 m~ 3 and Ks ~ 10~ 34 kg~ x ms 2 . With zero 
cosmological constant the black hole adiabatic compressibility at zero tem- 
perature is given by fl90l) with P = 0, 



«S| T =P=Q = — = n.8 » ( 94 ) 



25 _ 4ttM 2 G' 3 
"9~ ~ 9 c 8 

where the relevant factors of c and G are included, and the entropy has been 
set to the extremum value of 2irM. Putting in the numbers 

k s | t=0 = 2.6 x 10- 38 (J^-) kg-'ms 2 , (95) 



which is four orders of magnitude less than that of a solar-mass neutron star. 
We conclude that the zero temperature black hole equation of state, although 
"softer" than that of a non-rotating black hole, is still very much stiffer than 
that of a neutron star. 

The "softest" compressibility for a neutral black hole however is the 
isothermal compressibility: for an extremal black hole 

_ 25(11 + 80PS+ 128(PS) 2 ) 22 S 

KtIt=0 ~ (1 + APS){3 + A8PS + 128{PS) 2 ) "3~' (%) 
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some 33 times larger than Ks\ T=p=0 in f[54"j) . but still much larger than de- 
generate matter in a solar-mass neutron star. 

7 Discussion 

The obvious open question arising from the ideas presented here is: what 
about A > 0? 

The analysis of critical behaviour in §5 is only valid for A < 0, this 
critical point lies deep in the region P > and does not appear to be of 
any relevance to astrophysical situations. It is certainly of interest in the 
AdS-CFT correspondence [IB] but the particular analysis of §5, being in 
1 + 3-dimensions could only be relevant to 2 + 1-dimensional conformal field 
theory, which is of course of interest in its own right [H]. One could perform 
a similar analysis for 4 + 1-dimensional, or yet higher dimensional black holes, 
to try and gain insight into higher dimensional conformal field theory, and 
indeed this seems to have been the motivation in [421 14"3"1 [T4] . but the analysis 
of angular momentum in higher dimension is more difficult and will be left 
for future work. 

The thermodynamics of black holes in de Sitter space-time is a notoriously 
difficult problem [25[ [281 [291 [501 EH E21 [531 [51] as there are two event horizons 
and no "asymptotically de Sitter" region inside the cosmological horizon. 
Even with no black hole, a naive interpretation of the cosmological horizon 
implies that the transition from A = to any infinitesimally small A > 
appears to involve a discontinuous jump from zero to infinite entropy, at least 
if one associates the usual Hawking-Bekenstein entropy with the cosmological 
horizon when A > 0. 

Nevertheless it is argued in [5U] that a consistent strategy is to fix the 
relevant components of the metric at the cosmological horizon, rather than 
at spacial infinity as would be done in asymptotically flat or AdS space-time. 
When that is done the same expression for the ADM mass (121"]) is obtained, 
but with L 2 — > —L 2 , so H > 1 while the angular momentum is still given by 
J = aM. In this picture, all of the formulae in §3 are applicable for positive 
A and negative P, provided P is not too negative. If A is too large the black 
hole horizon and the cosmological horizon coincide and demanding that this 
does not happen puts a lower bound on P, for any fixed S, J and Q: with 
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Q = 0, for example, this requirement constrains P to 

p > ^f~* s . m 

Provided P lies above this lower bound we can analytically continue ( 1261) to 
negative P, with the understanding that S is the entropy of the black hole 
event horizon only and does not include any contribution from the cosmo- 
logical horizon. 

Of course P < is thermodynamically unstable, but it can be argued, 
in some circumstances at least, that positive pressures can be analytically 
continued to negative pressures [201 US], an d in a cosmological context there 
can now be little doubt that P < 0. Adopting the strategy of [50J the 
maximal efficiency of a rotating black hole in de Sitter space will be less than 
in the A = case, based on simply changing the sign of A in §4, and the zero 
charge efficiency vanishes when the black hole horizon and the cosmological 
horizon coincide at PS = — |. Any such deviation from the A = case 
will however be completely negligible for astrophysical black holes around 
one solar mass and the observed value of A, but it could be more significant 
during periods of inflation when A was larger. 

It has been suggested that primordial black-holes may have formed in 
the early Universe [S3] and, if this is the case and if they formed in sufficient 
numbers at any stage, then one should model the primordial gas as containing 
a distribution of highly incompressible black holes, like beads in a gas. These 
would certainly be expected to affect the overall compressibility of the gas 
as well as the speed of sound through the gas. In a radiation dominated 
Universe, ignoring the matter density, the speed of sound in the photon gas 
would be given by 

de 



c ~ 2 = 

7 op 



3c" 2 , (98) 



s 



where e is the energy density (essentially since the equation of state is P = ~e) 
so c 7 = 0.577 c. Since the speed of sound associated with the embedded black 
hole "beads" is Cs > y0^9 c = 0.9487 c the presence of a significant density of 
primordial black holes would expected to affect speed of sound in the photon 
gas and thus affect the dynamics. 

The square of the speed of sound for an extremal electrically neutral black 
hole is plotted in figure 5, for PS > —1/8. For comparison the asymptotic 
value (c| = 1/2 for PS — > oo) and the speed of sound in a thermal gas of 
photons (q? = 1/3) are also shown. 
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Figure 5: The speed of sound for an electrically neutral, extremal, black hole 
(with c = 1). 

8 Conclusions 

In conclusion there are strong reasons to believe that the cosmological con- 
stant should be included in the laws of black hole thermodynamics as a 
thermodynamic variable, proportional to the pressure of ordinary thermo- 
dynamics. The conjugate variable is a thermodynamic volume (J7|) and the 
complete first law of black hole thermodynamics is now (J5j), 

dU = TdS + VLdJ + $rfQ - PdV. (99) 

With this interpretation the ADM mass of the black hole is identified 
with the enthalpy 

M = H(S,P,J,Q) = U(S,V,J,Q) + PV (100) 

rather than the internal energy, U, of the system. 

The inclusion of this extra term increases the maximal efficiency of a 
Penrose process: for a neutral black hole in asymptotically anti-de Sitter 
space the maximal efficiency is increased from 0.2929 in asymptotically flat 
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space to 0.5184 in the asymptotically AdS case. For a charged black hole the 
efficiency can be as high as 75%. A positive cosmological constant is expected 
to reduce the efficiency of a Penrose process below the asymptotically flat 
space value. 

This point of view makes the relation between asymptotically AdS black 
holes and the Van der Waals gas, first found in [121113], even closer as there is 
now a critical volume associated with the critical point. The thermodynamic 
volume then plays the role of an order parameter for this phase transition 
and the critical exponents take the mean field values, 

a = 0, = 1 7 = 1, 5 = 3. (101) 

While there is no second order phase transition for a black hole in de 
Sitter space, there are other possible physical effects of including the PdV 
term in the first law. The adiabatic compressibility can be calculated ( 1901) 
and the speed of sound for such a black hole fl92|) is greater even than that 
of a photon gas and approaches c when PS = —1/8. 

Despite much progress the thermodynamics of black holes in de Sitter 
space-time is still very poorly understood and no doubt much still remains 
to be discovered. 
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